August 1997 


N A I R ET A L . 


UA&I- 5 L Q 


1753 


. * _ y ^7 


207278 


Numerical Simulation of the 9-10 June 1972 Black Hills Storm Using CSU RAMS 

U. S. Nair and Mark R. Hjelmfelt 

Institute of Atmospheric Sciences , South Dakota School of Mines and Technology . Rapid City, South Dakota 


Roger A. Pielke Sr. 

Department of Atmospheric Science, Colorado State University, Fort Collins, Colorado 
(Manuscript received 28 March 1996, in final form 4 December 1996) 

ABSTRACT 

Strong easterly flow of low-level moist air over the eastern slopes of the Black Hills on 9-10 June 1972 
generated a storm system that produced a flash flood, devastating the area. Based on observations from this 
storm event, and also from the similar Big Thompson 1976 storm event, conceptual models have been developed 
to explain the unusually high precipitation efficiency. In this study, the Black Hills storm is simulated using the 
Colorado State University Regional Atmospheric Modeling System. Simulations with homogeneous and inhom- 
ogeneous initializations and different grid structures are presented. The conceptual models of storm structure 
proposed by previous studies are examined in light of the present simulations. 

Both homogeneous and inhomogeneous initialization results capture the intense nature of the storm, but the 
inhomogeneous simulation produced a precipitation pattern closer to the observed pattern. The simulations point 
to stationary tilted updrafts, with precipitation falling out to the rear as the preferred storm structure. Experiments 
with different grid structures point to the importance of removing the lateral boundaries far from the region of 
activity. Overall, simulation performance in capturing the observed behavior of the storm system was enhanced 
by use of inhomogeneous initialization. 


1. Introduction 

The 9-10 June 1972 Black Hills storm resulted in the 
most devastating flash flood in the recorded history of 
South Dakota. The storm was characterized by convec- 
tive cells of high precipitation efficiency characteristic 
of tropical precipitation systems rather than that of a 
typical midlatitude convective system. The intense pre- 
cipitation from the orographically forced storm system 
over the basins of Rapid, Boxelder, Spring, and Battle 
Creeks (Schwarz et al. 1975) led to severe flooding, 
resulting in 237 deaths and property damage in excess 
of $160 million (see Fig. 1). 

Large-scale features associated with the Black Hills 
storm and the Big Thompson (Colorado) storm of 1 
August 1976, which was similar in nature, are sum- 
marized by Maddox et al. (1978). Conceptual models 
that explain high precipitation efficiency were suggested 
for the Black Hills storm by St. Amand et al. (1972) 
and for the Big Thompson storm by Caracena et al. 
(1979). Dennis et al. (1973) discuss the motion of con- 
vective cells observed for the Black Hills storm. Nu- 
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merical simulations of the Black Hills storm were pre- 
viously made using a two-dimensional, time-dependent 
cloud model (Kopp and Orville 1973). More recently, 
mesoscale simulations have been made by Chang and 
Perkey (1995). Yoshizaki and Ogura (1988) investigated 
the Big Thompson storm using two- and three-dimen- 
sional cloud models. This study focuses on testing the 
ability of the Colorado State University Regional At- 
mospheric Modeling System (CSU RAMS) to capture 
the essential nature and fury of the storm and to examine 
the various hypotheses for structure and motion of the 
storm in light of the simulation results. 

Section 2 describes the observations of the storm 
event. The hypotheses about the storm structure for the 
Black Hills storm (St. Amand et al. 1972) and the Big 
Thompson storm (Caracena et al. 1979) will be exam- 
ined in section 3. A brief description of the model in- 
itialization is provided in section 4. The results from 
the model simulations are described in sections 5, 6, 
and 7, and the study is summarized in section 8. 

2. Observations 

The 0000 UTC surface chart for 10 June 1972 (Fig. 
2) shows a polar front to the southeast of the Black Hills 
region. To the rear, east of the front, a band of moist 
air with dewpoints ranging from 18° to 21°C can be 
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Fig. 1. (a) Location of area affected by the storm; (b) topography for the Black Hills region. 
Also shown are the major creeks involved in the flooding, and the location of Rapid City is 
indicated by the cross. 


observed. The dominant feature of the wind fields in 
the surface chart is strong easterlies of 10-15 m s -1 to 
the east of the front. 

A strong negatively tilted ridge over South Dakota is 
the main feature observed on the 0000 UTC 500-mb 
chart (Fig. 3). A weak trough can be seen approaching 
the Black Hills region from the southwest. The 0000 
UTC 10 June 1972 Rapid City sounding (Fig. 4) shows 
strong east-southeasterly flow through a deep layer 


above the surface. The sounding exhibits high instability 
with a lifted index of -6. The wind veers to light west- 
erlies in the upper levels. The water vapor mixing ratio 
in the lower layer is of the order of 14 g kg '. The layer 
immediately above the surface is conditionally unstable. 
If the surface air is forced upward, it could become 
saturated and positively buoyant above 750 mb. 

The 1200 UTC 9 June 1972 sounding for Huron (fig- 
ure not shown), which is characteristic of the air mass 
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Fig. 3. The 500-mb map for 0000 UTC 10 June 1972 (Schwarz et ah 1975). Bold solid lines are height 
contours in dekameters and the bold dashed lines are isotherms in degrees Celsius. 
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Fig. 4. Rapid City sounding for 0000 UTC 10 June 1972 (Dennis 
et aL 1973). 


to the rear of the front, showed strong temperature in- 
versions. These inversions prevented the release of in- 
stability until the air mass experienced lifting over the 
slopes of the Black Hills. The shortwave trough helped 
in the release of instability and also, toward the end of 
the storm event, to flush the convective cells onto the 
eastern plains. 

The track of radar echoes before 0000 UTC (Fig. 5), 
as observed by a research X-band radar operated by the 
Institute of Atmospheric Sciences (Dennis et al. 1973), 
shows cells forming to the south of the Black Hills and 
moving northwest with time. However, the long-term 
trend of the radar echoes (Fig. 6) was to remain on the 
eastern slopes of the Black Hills. Figure 6 is a composite 
drawing by Schwartz et al. (1975) of significant radar 
echoes for the period 0000-0700 UTC, based on ob- 
servations by the National Weather Service (NWS) ra- 
dar at Huron, South Dakota (300 km east of Rapid City), 
and the Air Force radar at Ellsworth Air Force Base (15 
km east of Rapid City). 

The total accumulated rainfall from the storm event 
is shown in Fig. 7. The rainfall pattern is aligned in a 
north-south direction to the west of Rapid City. The 
maximum observed point accumulation from the storm 
event was 380 mm. The area of heavy rainfall was re- 
stricted to the eastern slopes of the Black Hills, between 
Rapid City and the central ridge of the hills. 


3. Hypothesis on storm motion and intensity 

Based on these observations, models have been pro- 
posed as to the structure of the storm and also the motion 
of convective cells. Based on radar observations prior 
to 0000 UTC 10 June, Dennis et al. (1973) suggest a 
long succession of convective cells moving from the 
southeast and suddenly intensifying to drop heavy show- 
ers on the eastern Black Hills and rapidly diminishing 
as they continue northwest or north-northwest. Dennis 
et al. (1973) also note the pronounced absence of down- 
drafts during the storm event. 

Dennis et al. (1973) estimate the precipitation effi- 
ciency of the storm system between 80% and 90%. To 
account for the high precipitation efficiency of the sys- 
tem, St. Amand et al. (1972) invoke recycling of water 
substance between mature and newly forming cells. The 
St. Amand et al. (1972) model (Fig. 8) has a mature 
cloud mass over the Black Hills due to lifting of moist 
easterly flow that slopes strongly to the east above 450 
mb. Precipitation from the upper sheared region of ma- 
ture cloud mass is introduced into newly formed cells 
downstream. This recycling of precipitation enhances 
the precipitation efficiency of the storm system. 

Based on observations from the Big Thompson storm 
of 1976, Caracena et al. (1979) proposed a conceptual 
model of the storm structure that would explain the high 
precipitation efficiency and quasi-steady state. The pre- 
cipitation efficiency for the Big Thompson storm was 
estimated to be around 85%. According to Caracena et 
al. (1979), for orographic convective storms in reverse 
shear, the strong boundary layer shear underlying weak 
upper-level winds roots the lower portion of cloud on 
the mountain. The vertical transport of easterly mo- 
mentum causes the updraft structure to tilt westward 
with height. The resulting structure (Fig. 9) allows the 
precipitation to fall to the rear of the updraft. This en- 
ables the updraft to be relatively free of precipitation 
loading and thus exist in a quasi-steady state. 

Numerical simulation of the Black Hills storm by 
Kopp and Orville (1973) using a two-dimensional time- 
dependent model showed a tendency for the storm to 
drift over the hills and out of the domain (westward 
drift). They increased the strength of the upper-level 
winds, which countered this tendency of westward drift. 
The numerical simulation of the Big Thompson flood 
by Yoshizaki and Ogura (1988) showed storm structures 
similar to that proposed by Caracena et al. (1979). Their 
experiments using homogeneous initialization of mois- 
ture resulted in a region of intense precipitation located 
at the peak of the terrain, contrary to observations. They 
also noted the formation of a cold pool underneath the 
convective cells. 

4. Model description and initialization 

The CSU RAMS is a mesoscale cloud model that has 
both hydrostatic and nonhydrostatic formulations. In- 
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teractive nested grid structure is available to simulate 
phenomena that require finer grid spacing. In these sim- 
ulations, the model utilizes a bulk water parameteriza- 
tion for representing cloud microphysics. A detailed de- 
scription of RAMS can be found in Pielke et al. (1992). 

For the present study, the nonhydrostatic formulation 
was utilized. Four experiments were performed, three 
with horizontally homogeneous initialization and one 
with inhomogeneous initialization. For both modes of 
initialization, 0000 UTC data were used. The horizon- 
tally homogeneous grid used in the first experiment (re- 
ferred to from here on as HH1) was initialized using 
the 0000 UTC Rapid City sounding. The model domain 
(Fig. 7), spanned most of the region of the Black Hills 
where the greatest amount of precipitation was observed 
(90 km X 150 km region centered at - 103.5744°). The 
grid spacing utilized was 1 km in the horizontal. In the 
vertical, the grid spacing was 100 m near the surface, 
with a grid stretch ratio of 1.15 and maximum vertical 
spacing of 1 km above 7 km. The time step used for 
the horizontally homogeneous simulation was 10 s. 


The inhomogeneous initialization experiment (re- 
ferred to from here on as IH1) utilized a two-way in- 
teracting nested-grid configuration. Three grids were 
used to nest down to the region of interest in the Black 
Hills. The outermost, coarse-grid spacing (grid 1) of 18 
km covered most of South Dakota and parts of bordering 
states (North Dakota, Wyoming, Nebraska, Montana). 
Nested within this grid was a finer grid (grid 2) of 6-km 
grid interval that covered most of the Black Hills region. 
The finest grid (grid 3), nested in the second grid, had 
a grid interval of 2 km and spanned the region of the 
Black Hills where most of the precipitation from the 
storm fell (Fig. 7). The vertical grid spacing for all three 
grids was 100 m near the surface, with a grid stretch 
ratio of 1.2 and maximum grid spacing of 750 m above 
4 km. Variable initialization utilized soundings and sur- 
face observations at 0000 UTC from South Dakota and 
neighboring states. The time steps for grids 1, 2, and 3 
were 30, 15, and 5 s, respectively. The model topog- 
raphy was derived from the 30" U.S. Navy terrain da- 
taset. Figure 10 shows the topography for the horizon- 
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Fig. 6. Composite of radar echoes for the period of 0000-0700 UTC 10 June 1972 (Schwarz et al. 1975). 


tally homogeneous (HH1) grid and grid 3 from the in- 
homogeneous initialization grid. The boundary condi- 
tions used for both the HH1 and IH1 runs are a rigid-top 
boundary and the Klemp-Wilhelmson (1978) radiation 
boundary condition for lateral boundaries. For the in- 
homogeneous initialization simulation, the nudging op- 
tion that allows nudging of boundaries toward the next 
12-h observation was also activated. In the outer two 
grids for the inhomogeneous initialization run, the mi- 
crophysical parameterization was not activated and wa- 
ter vapor was treated as a tracer. 

Two additional experiments with horizontal-homo- 
geneous initialization were made using the following 
grid structures: 1) single grid with the same vertical and 
horizontal structure as that of the innermost grid used 
in IH1; 2) nested grid structure, the same as that used 
in IH1. These experiments will be referred to from here 
on as HH2 and HH3, respectively. Except for the vertical 
and horizontal grid spacing, all the other model settings 
are the same for HH1 and HH2. Experiments HH3 and 
IH1 differ only in that for HH3 all the grids are ini- 
tialized using the 0000 UTC Rapid City sounding and 
the nudging option is not used. 

5. HH1 experiment results 

The model was integrated for 6 h, starting from initial 
conditions at 0000 UTC. Clouds started developing after 
15 min of integration as the moist conditionally unstable 
air was forced up the slopes by low-level easterly flow. 


In early stages of the simulation, clouds form in bands 
oriented in a north-south direction in the center of the 
grid. At later times, clouds drift to the left side of the 
model grid, that is, toward the west. An hour into the 
simulation, cloud fields start exhibiting a pattern of new 
cells developing in the southwestern part of model grid 
with individual cells propagating in a northwest direc- 
tion. The cells align themselves parallel to the flow. 
Figure 1 1 shows clouds at 15-min intervals starting from 
an hour of model integration, illustrating this pattern. 
This pattern of movement persisted throughout the 6 h 
of simulation and is reminiscent of radar observations 
(Fig. 5) reported by Dennis et al. (1973) discussed in 
section 3. Precipitation is delayed until the cells have 
reached the higher elevations of the central hills. 

The simulated storm system exhibited a precipitation 
pattern different from that observed. The observed pre- 
cipitation pattern over the region covered by the grid 
(Fig. 12b) was oriented north-south approximately sym- 
metric about the center of the grid, while the simulated 
precipitation pattern is oriented with major axis south- 
east-northwest (Fig. 12a) and lies in the western region 
of the grid. The maximum accumulated precipitation at 
the end of 6 h was 275 mm, compared to an observed 
maximum of 380 mm. During the simulation, some of 
the cells drifted westward and moved out of the model 
domain. This tendency of westward drift was also noted 
in a modeling study by Kopp and Orville (1973). The 
region of intense precipitation for the simulation occurs 
at the high terrain region, in contradiction to the ob- 
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Fig. 7. Total accumulated precipitation for 9-10 June 1972 for a period of 16 h. Contour interval 
is 50 mm. Also shown are the locations of the grids used in the experiments. Rectangle labeled 
A is the location of the grid used in homogeneous initialization (HH1), and rectangle labeled B 
is location of grid 3 in inhomogeneous simulation (adapted from Schwarz et al. 1975). 
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Fig. 8. Conceptual model for storm structure proposed by St. Amand et al. (1972). 
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Fig. 9. Conceptual model for the storm structure proposed by Caracena et al. (1979). 


servation. Similar results were obtained in the modeling nificant strong downdrafts or cold outflows, particularly 

study of the Big Thompson flood by Yoshizaki and Ogu- in regions of precipitation shafts. This is consistent with 

ra (1988). One of their experiments that had horizontally the observations by Dennis et al. (1973) but is in contrast 

homogeneous initial distribution of moisture resulted in to the cold air pools beneath the storms reported by 

heavy precipitation occurring at the mountain peak, in- Yoshizaki and Ogura (1988) in their modeling study of 

consistent with the observations. the Big Thompson storm. The lack of downdrafts has 

Vertical cross sections through the convective cells been attributed to moist midlevel air and low cloud bas- 

(Fig. 13) show an updraft with very little tilt with height, es, which restrict evaporative cooling. Model results in- 

Another interesting feature is the absence of any sig- dicate graupel melting to form rain was an important 
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Fig. 10. The topography for (a) HH 1 and (b) grid 3 of IH1. Contours in intervals ot 100 m 
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Fig. 1 1. Cloud field at 3.7 km for the homogeneous initialization at (a) 60, (b) 75, and (c) 90 min into the 
model integration. The contour interval is 0.75 g kg 
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Fig. 12. (a) The 6-h accumulated precipitation for homogeneous model initialization (HHI) in millimeters; (b) total accumulated precipi- 
tation from the storm event tor the same area in millimeters; (adapted from Schwarz et al. 1975). Contour interval is 50 mm. 


microphysical process, but it appears that cooling 
achieved from the melting process was not strong 
enough to initiate strong downdrafts. 

6. Inhomogeneous initialization (IH1) experiment 
results 

The inhomogeneous simulation was integrated for 6 
h starting from the initial conditions. Clouds started 
forming near 30 min of simulated time, about 15 min 
later than in the HHI simulation. The storm structure 
and motion were significantly different than for the hor- 
izontally homogeneous simulation. The storms formed 
in the center of the model grid (grid 3) in a band aligned 
in a north-south direction and tended to be stationary. 
This is consistent with the long-term trend of radar ech- 
oes (Fig. 6) discussed in section 2. 

The vertical section through the cells (Fig. 14) shows 
a pattern very similar to the conceptual model proposed 
by Caracena et al. (1979). A very pronounced tilt of the 
updraft in the east-west direction with height can be 
seen. Also, the rainfield overlaid on the vertical velocity 
field shows the precipitation falling to the rear of the 
updraft. The IH1 simulation lacked any strong down- 


draft features, but compared to the HHI simulations it 
had some regions of weak downward motion. 

The accumulated precipitation field for the IH1 sim- 
ulation compares more favorably with the observed pre- 
cipitation pattern than did the HHI simulation. Most of 
the precipitation from the simulation is in the eastern 
half of the domain (Fig. 15a), as is the observed pre- 
cipitation pattern (Fig. 15b). Also, the simulation suc- 
ceeds in approximating the locations of “hot spots,” or 
areas of most intense precipitation. The maximum 6-h 
accumulated precipitation for the IH1 is 220 mm. 


7. HH2 and HH3 experiment results 

The HH2 and HH3 experiments were conducted to 
investigate substantial differences between the HHI and 
IH1 experiments. In the HH2 simulation (single grid but 
with a larger grid spacing, similar to that in grid 3 for 
the nested cases), the cloud patterns and motion were 
similar to that of HH 1 . The general pattern of 6-h ac- 
cumulated precipitation field (Fig. 16a) was similar to 
that obtained for the HHI experiment. The magnitude, 
however, was very different. The maximum value of the 
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Fig. 13. East-west cross sections at 55 km from southern boundary of (a) cloud held, contour interval of 
0.8 g k» (b) vertical velocity, rain mixing ratio (contours in white) and wind vectors for HH1 120 min 
into the simulation. The rain mixing ratio contour interval is 0.3 g kg 1 and the contour interval for vertical 
velocity held is 3 m s '. 


6-h accumulated precipitation for HH2 was 106 mm, 
compared to 275 mm obtained for the HH1 experiment. 

The HH3 experiment utilized the nested grid of IH1 
but with a homogeneous initialization. The nature of 
cloud patterns for the HH3 experiment was very similar 
to that of IH 1 . The 6-h accumulated precipitation pattern 
(Fig. 16b) was oriented in a north-south direction, sim- 
ilar to the IH1 experiment results. Experiments IH1 and 
HH3 have different values of maximum accumulated 
precipitation. The HH3 experiment has a maximum val- 
ue of 305 mm, compared to 220 mm obtained for the 
IH1 experiment. The HH3 experiment, however, con- 
centrates the precipitation to a single maximum, com- 
pared to four local maxima for the IH1 experiment. The 
precipitation field of the IH1 experiment compares bet- 
ter with the observed precipitation pattern in that re- 
spect. 


8. Conclusions 

Based on comparison of the observed and simulated 
precipitation patterns, it could be concluded that the 
inhomogeneous simulation (IHl) performed better than 
the horizontally homogeneous simulations. Both types 
of simulations rule out the St. Amand et al. (1972) con- 
ceptual model for the structure of the Black Hills 9 June 
1972 storm. The IHl simulation results strongly support 
the conceptual model suggested by Caracena et al. 
(1979) for the Big Thompson flood storm structure as 
also being applicable to the Black Hills storm. Homo- 
geneous simulation results are consistent with early ob- 
servations (before 0000 UTC 10 June 1972) of the cell 
motions reported by Dennis et al. (1973), but the long- 
term trend of radar echoes is more consistent with the 
IHl simulation. 
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Fig. 14. East-west cross sections at 90 km from the southern boundary of (a) cloud field, contour interval 
of 0.8 g kg (b) vertical velocity, rain mixing ratio (contours in white) and wind vectors for IHI 120 min 
into the simulation. The rain mixing ratio contour interval is 1.5 g kg ' and the contour interval for vertical 
velocity held is 6 m s 


In our preliminary study (Nair et al. 1995), we at- 
tributed the differences between the simulations HH1 
and IHI to differences in the spatial distribution of mois- 
ture between the two initializations. It was proposed that 
IHI succeeded in including the air mass to the east of 
the Black Hills, thus incorporating the characteristics of 
air mass in which the storm actually developed after 
0000 UTC 10 June 1972. The experiments HH2 and 
HH3 were conducted to investigate the effect of bound- 
ary conditions and grid spacing on the solution. The 
HH2 grid, which has the same structure as that of the 
innermost grid of the IHI nested grid, showed similar 
behavior to the HH1 experiment. The larger grid spacing 
used in HH2 resulted in lower values of 6-h accumulated 
precipitation. Thus, the difference in grid spacing cannot 
explain the differences in the nature of evolution of 
storms in the HH1 and IHI experiments. The HH3 ex- 


periment, which has the same nested grid structure and 
boundary conditions as the IHI experiment (except for 
nudging), resulted in storm evolution and a precipitation 
pattern similar to that obtained in IH 1 . This would imply 
that the differences between HH1 and IHI are primarily 
due to the differences in boundary conditions. The nest- 
ed grid structure, which moves the artificial lateral 
boundary away from the region of interest, provides 
results in better agreement with the observations. The 
pattern of the precipitation field from the IHI experi- 
ment provides the best match with the observed pre- 
cipitation field. However, the maximum values of ac- 
cumulated precipitation for the HH3 simulation are clos- 
er to the observations. Since the spatial distribution of 
precipitation is important for flood events, and the mag- 
nitude of accumulated precipitation from the IHI sim- 
ulation is high enough to cause flooding (Kenner et al. 
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1996), we conclude that the IH1 simulation compares 
most favorably with observations. Inhomogeneous in- 
itialization with nested grid structure is preferred for 
good predictions in situations similar to the present 
study. 
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